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ABSTRACT 

Planetesimal accretion during planet formation is usually treated as collision- 
less. Such accretion from a uniform and dynamically cold disk predicts proto- 
planets with slow retrograde rotation. However, if the building blocks of proto- 
planets, planetesimals, are small, of order of a meter in size, then they are likely 
to collide within the protoplanet's sphere of gravitational influence, creating a 
prograde accretion disk around the protoplanet. The accretion of such a disk 
results in the formation of protoplanets spinning in the prograde sense with the 
maximal spin rate allowed before centrifugal forces break them apart. As a result 
of semi-collisional accretion, the final spin of a planet after giant impacts is not 
completely random but is biased toward prograde rotation. The eventual accre- 
tion of the remaining planetesimals in the post giant-impact phase might again 
be in the semi-collisional regime and delivers a significant amount of additional 
prograde angular momentum to the terrestrial planets. We suggest that in our 
Solar System, semi-collisional accretion gave rise to the preference for prograde 
rotation observed in the terrestrial planets and perhaps the largest asteroids. 

Subject headings: planets and satellites: formation — solar system: formation 

1. INTRODUCTION 

Protoplanets form by the accretion of planetesimals. When planetesimals are accreted, 
they deliver rotational angular momentum due to their relative motion with respect to the 
protoplanet. This accretion is usually treated as collisionless, assuming that collisions among 
planetesimals can be neglected while they are within the Hill sphere of the protoplanet. In 
collisionless accretion, the angular momentum accreted from a uniform and dynamically cold 
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disk of planetesimals results in slow retrograde rotation (Lissauer & Kary 1991; Dones & 
Tremaine 1993a). Lissauer et al. (1997) have shown that rapid prograde rotation can only 
be achieved if disk density profiles are imposed such that the surface mass density near the 
outer edges of a protoplanet's feeding zone is significantly greater than that in the rest of the 
accretion zone. This suggests that protoplanets do not possess any significant spin due to 
collisionless planetesimal accretion. The final stage of terrestrial planet formation consists 
of collision and accretion events of a few dozen protoplanets of about 0.05 M e (Agnor et al. 
1999; Chambers 2001; Goldreich et al. 2004b). These giant impacts deliver spin angular 
momentum to the final planet. If giant impacts are solely responsible for the final spin 
properties of terrestrial planets then terrestial planets should display random obliquities 
(the angle between the orbital and rotational angular momentum) and exercise prograde 
and retrograde rotation with equal likelihood. 

In this paper, we investigate the possibility of semi-collisional or collisional planetesimal 
accretion and the effect it would have on planetary spins. In §2 we first determine the range 
of planetesimal sizes for which semi-collisional or collisional accretion applies and derive 
the consequences of semi-collisional and collisional accretion for the spin of protoplanets. 
The spin of terrestial planets due to giant impacts of protoplanets is calculates in §3 and 
compared with the semi-collisional contribution. Post giant-impact accretion is discussed in 
§4. Comparison with the Solar system and conclusions follow in §5. 

2. SEMI-COLLISIONAL AND COLLISIONAL ACCRETION 

The Hill radius denotes the distance from the protoplanet at which the tidal force due 
to the Sun and the gravitational force due to the protoplanet both acting on a planetesimal 
are in equilibrium. It is given by 



where a is the semi-major axis of the protoplanet and m its mass. When two planetesi- 
mals collide with each other while passing through the Hill sphere of the protoplanet, one 
or both of them become bound to the protoplanet. Further collisions among the bound 
particles damp their random motions, leading to the formation of an accretion disk around 
the protoplanet (Sari & Goldreich 2006) (see figure 1). Inelastic planetesimal collisions and 
subsequent capture by the planet's gravitational field has been proposed in order to form 
circumplanetary disks from which regular satellites could form (e.g. Safronov et al. 1986; 
Estrada & Mosqueira 2006). Here we are exploring the possibility that the growth of pro- 
toplanets is dominated by the accretion of such a planetesimal disk. The details of this 
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accretion process, such as what fraction of bound particles will be accreted by the growing 
protoplanet, are uncertain. Perturbations from nearby protoplanets and moons or gas, if still 
present at the time protoplanets form, may facilitate the dissipation of the planetesimals' 
angular momentum allowing efficient accretion onto the protoplanet. 

2.1. Planetesimal Sizes 

For dynamically cold planetesimal disks, the ratio between the rate of planetesimal 
collisions within the Hill sphere to the rate of direct collisions onto the protoplanet is r 9 « -1 / 2 , 
where r g is the optical depth within the disk plane over a distance of Rh and a = r / Rh where 
r is the protoplanet's radius. If r g a~ x l 2 > 1, the accretion may be dominated by binding 
planetesimals into an accretion disk rather than direct impacts onto the protoplanet; we 
call this semi-collisional accretion. Collisional accretion takes over for r g > 1. For inelastic 
planetesimals with velocity u < vh = QRh, the optical depth in the disk is given by 
r g ~ 3<jvh I 'sp s u where Q, s, p s and a are the protoplanet's Keplerian angular velocity 
around the Sun, the typical planetesimal radius, material density and overall mass surface 
density respectively. The random velocities of the planetesimals are damped by mutual 
collisions and stirred by gravitational interactions with the protoplanets. When these two 
processes are in equilibrium, we have 

— ~ a , tor u < vh (2) 

vh cr r 

where £ is the mass surface density of the protoplanets (Goldreich et al. 2004b). Most of 
the planetesimal accretion occurs when S/u ~ 1. The condition for semi-collisional accre- 
tion {jgor^l 2 > 1) together with equation 2, defines an upper limit to the planetesimal size 
for which semi-collisional accretion holds. Using the minimum mass solar nebula (Hayashi 
1981) surface density of ~ 8 g/cm 2 at 1 AU, p s ~ 3 g/cm 3 and an isolation mass ~ 0.05 M® 
(Weidenschilling et al. 1997), we find 1 that s < 9 m. A lower limit to the planetesimal size 
is given by the velocity dispersion for which the disk becomes locally unstable to gravita- 
tional collapse. This velocity is ~ 10 cm/s at 1 AU, corresponding to a minimum size for 
planetesimals of ~ 6 cm. Therefore semi-collisional or collisional accretion applies as long as 
6 cm < s < 9 m. A fragmentation cascade produced by destructive planetesimal collisions 
leads to the formation of ever smaller planetesimals (Goldreich et al. 2004b). In fact, grav- 
itational instabilities in the disk may be responsible for the lower limit on the planetesimal 



1 All estimates above assumed u < o^^vu- However, for large enough planetesimals we have oi^^vh < 
u < vh- Taking this into account results in a slightly higher upper limit of 17 m for s. For simplicity and 
given the order of magnitude nature of this calculation, we ignore this complication. 
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size, in which case s ~ 6 cm. Possible gaseous remnants of the solar nebula may lower the 
velocity dispersion preventing fragmentation down to the stability limit. Though this is an 
uncertainty during protoplanet formation it is unlikely that significant amounts of gas pre- 
vailed after giant impacts. Further, the low bulk density (~ 0.6 g/cm 3 ) of comets (A'Hearn 
et al. 2005; Davidsson & Gutierrez 2006) seems to suggest gentle accretion of small bodies 
and therefore supports the idea of semi-collisional or collisional accretion. 



2.2. Spin of Protoplanets due to Planetesimal Accretion 

We assume that the orbits of the planetesimals and the protoplanets are circular and 
co-planar. The interaction between the planetesimals and the protoplanet can be described 
by Hill's equations (Hill 1878; Goldreich & Tremaine 1980; Petit & Henon 1986). In our 
coordinates the position of the planetesimal is given with respect to the protoplanet. The 
x-axis points radially outwards and the y-axis in the prograde direction. The equations of 
motion are given by 

Gun 

x-2Qy-3Q 2 x = - {x2+y2)3/2 x (3) 

„^ . Gra , 4 . 

y + 2Qx = — y. 4 

(ar + y 1 ) 6 ! 1 

We solve these equations numerically and sum the specific angular momentum of all plan- 
etesimals that pass within some effective accretion radius R aC c- In collisionless accretion the 
protoplanet accretes at its actual radius so that R acc = r; in semi-collisional or collisional 
accretion an accretion disk forms and the protoplanet effectively accretes at its gravitational 
radius such that R acc ~ Rh- Figure 2 shows that protoplanets acquire a retrograde spin 
for R acc < 0.2 R H and a prograde rotation for R acc ~ R H . The prograde rotation for 
Race » Rh can be understood by considering the angular momentum supplied by plan- 
etesimals solely due to the Keplerian shear of the disk. In this case the specific angular 
momentum acquired by the protoplanet is R^ cc fl/4 in the prograde sense (Lissauer & Kary 
1991; Dones & Tremaine 1993a). The actual angular momentum delivered to the planet is 
given by figure 2 for collisionless accretion only. In the semi-collisional and collisional cases, 
the disk must lose angular momentum before it can be accreted by the protoplanet. The 
accretion of such a disk results in the formation of protoplanets spinning in the prograde 
sense with the maximal spin rate allowed before centrifugal forces break them apart. 
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3. GIANT IMPACTS 

The final stage of terrestrial planet formation consists of collision and accretion events 
among the protoplanets. These giant impacts deliver spin angular momentum to the final 
planet. Provided that the random velocities of the protoplanets are sufficiently large, one 
can neglect the shear imposed by the differential rotation of the disk, so there is no preferred 
direction for giant impacts to occur. Giant impacts therefore deliver angular momentum 
in a random walk like fashion. Lissauer & Safronov (1991) and Dones & Tremaine (1993b) 
calculated the magnitude of the random component of the spin angular momentum due to a 
single giant impact and compared it with the observations. Here we determine the random 
and systematic spin angular momentum delivered to the final planet by N giant impacts 
using the following toy model. We start with N + 1 identical protoplanets all of mass m and 
radius r which are sequentially accreted one by one. After N such accretion events, we are 
left with a final planet of mass M = (N + l)m and radius R = (N + l) 1 / 3 r. We assume 
throughout that protoplanets are spherical with constant density p. 



3.1. Random Component of the Angular Momentum 

In the limit that the protoplanets' velocity dispersion is small compared to their impact 
velocity and assuming that protoplanets have no spin, the maximum angular momentum 
delivered by one impact is 

M T m 



Ima* = tt^— a/2G (Mr + m ) Wt + r) (5) 

where Mt is the mass and Rt the radius of the target. The root mean square (rms) angular 
momentum in the direction perpendicular to the plane of the Solar System (z-direction) con- 
tributed by a single impact can be obtained by averaging over all possible impact parameters 
and is given by l z rms = y/l/6lmax- Adding the contributions of each impact in quadrature, 
with Mt = nm and Rt = n 1 ^ 3 r for n — 1, 2, N, the final rms angular momentum in the 
^-direction after iV 3> 1 impacts is 

Kms = ^N'^ruMR 2 (6) 

where 



The precise number of giant impacts during the late stage of planet formation is uncertain. 
However, the final "isolation" mass for the minimum mass solar nebula at 1 AU is about 
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0.05 M e (Weidenschilling et al. 1997; Goldreich et al. 2004a). This suggests that about 20 
giant impacts have occurred in order to form an earth at 1 AU. For iV ~ 20, equation 6 
predicts a spin period of ~ 4 hours for the Earth. N-body simulations find a somewhat 
shorter spin period of ~ 1.8 hours for bodies more massive than 0.5 M e (Agnor et al. 1999). 
This rapid rotation originates from unphysical mergers between protoplanets encountering 
each other at more than the escape velocity. As expected, N-body simulations also show that 
final obliquities due to giant impacts with no initial spin are randomly distributed (Agnor 
et al. 1999; Chambers 2001). 



The final spin of a terrestrial planet after giant impacts is no longer random but contains 
a systematic component if each protoplanet possess a systematic spin due to semi-collisional 
planetesimal accretion. The systematic component of the angular momentum delivered by 
N impacts of maximally spinning protoplanets with prograde rotation is 



Comparing the random z-component of the angular momentum (equation 6) to the 
systematic one (equation 8), we find that they are similar in magnitude with the random 
component up to twice the systematic one for 1 < N < 60. The final distribution for the z- 
component of the angular momentum is obtained by combining the random and the ordered 
contributions. It is normally distributed with its mean given by equation 8 and its standard 
deviation given by equation 6. Since the mean is positive, corresponding to prograde rotation, 
we expect more prograde than retrograde spins in a given planetary system. We find about 
70% of all planets to be rotating in the prograde sense and only 30% in a retrograde manner 
for 10 < N < 60 giant impacts. 



3.2. Systematic Component of the Angular Momentum 



Spin — J^Spin 




(8) 



3.3. Comparison 



3.4. 



Uncertainties 



The following uncertainties could affect our estimates for prograde and retrograde rota- 
tion. We have assumed that the velocity dispersion of the protoplanets is small compared to 
the impact velocity. However, the velocity dispersion might be as large as the escape veloc- 



7 



ity from the protoplanet, in which case the random component of the angular momentum 
could increase up to ~ y/2. A higher fraction of planets with retrograde rotation would be 
produced if the mutual accretion of protoplanets were pairwise, such that all giant impacts 
were between equal-sized bodies, rather than one by one. Furthermore, the majority of the 
mass accreted is likely due to collisions close to head on, which deliver a smaller random 
component of angular momentum than grazing ones (Agnor & Asphaug 2004). On the other 
hand, grazing collisions could deliver spin angular momentum and little mass, allowing eas- 
ily twice the naive number of giant impacts. These uncertainties could be addressed using 
hydrodynamic simulations. 



The stirring force protoplanets exert on each other can be balanced by the force due 
to dynamical friction caused by the planetesimals as long as a > E, ensuring small random 
velocities of the protoplanets. However, as the protoplanets accrete more planetesimals, 
their surface density increases and dynamical friction becomes less and less effective until 
it is no longer able to balance the mutual stirring of the protoplanets. Orbit crossing and 
giant impacts set in when a ~ E (Goldreich et al. 2004a). Planetesimal accretion continues 
and additional 'new' planetesimals are produced as byproducts of giant impacts. The exact 
amount of smaller particles produced in a giant impact depends on the mass ratio of the two 
colliding protoplanets, their relative velocity and impact angle. For example, for collisions 
between like-sized protoplanets with an impact velocity of twice their escape velocity and an 
impact angle of 30° (where 0° corresponds to a head-on collision) about 10% of the total mass 
of the system escapes as smaller particles (Agnor & Asphaug 2004). Due to the production 
of 'new' planetesimals in giant impacts and the fact that giant impacts set in when a ~ E, 
large amounts of planetesimals are expected to still be present after the culmination of giant 
impacts. This is also required to relax the high eccentricities of planets expected after giant 
impacts (Goldreich et al. 2004a). N-body simulations predict eccentricities of ~ 0.1 for 
terrestrial planets after giant impacts (Chambers 2001). The eccentricity damping timescale 
tdamp due to dynamical friction caused by left-over planetesimals is given by 



This timescale should be shorter than the time required for the remaining planetesimals to 
be accreted onto the terrestial planets: 



4. 



ACCRETION AFTER GIANT IMPACTS 




(9) 




(10) 



8 



This yields 




2 



<7>E 



0.07S. 
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Therefore, more than 7% of the mass should still reside in planetesimals in order to damp the 
planets' eccentricities. The eventual accretion of the remaining planetesimals delivers addi- 
tional angular momentum to the planet. For sufficiently small planetesimals, this accretion 
would again be in the semi-collisional or collisional regime and hence deliver additional pro- 
grade angular momentum to the planet. The accretion of about 10% M e in semi-collisional 
manner would be sufficient to deliver an angular momentum equivalent to that of the Earth- 
Moon system. For Mars less than 3% of its mass would need to be accreted semi-collisionally 
to supply its current angular momentum, assuming that it had no previous spin. These small 
percentages indicate that semi-collisional or collisional accretion of only a small fraction of 
the planet's mass after giant impacts is sufficient to substantially alter planetary spins lead- 
ing again to favoritism of prograde rotation. Formation of gaps in the planetesimal disk after 
giant impacts may complicate this picture. 



We have shown here, that planetesimal accretion might be in the semi-collisional or col- 
lisional regime leading to the formation of a prograde accretion disk around the protoplanet. 
Such a disk gives rise to a maximally-spinning protoplanet with prograde rotation. The 
final spin of terrestrial planets is therefore no longer random but is biased toward prograde 
rotation. The dominance of prograde rotation might be increased further by the accretion of 
leftover planetesimals in the post-giant impact phase, provided that semi-collisional or colli- 
sional accretion still applies. Comparing our results with the spin properties of the terrestrial 
planets is somewhat difficult since the spins of Mercury and Venus have evolved consider- 
ably since their formation (Laskar & Robutel 1993; McCue & Dormand 1993), leaving only 
Earth and Mars whose spins have evolved to a much lesser degree. Earth and Mars both 
display prograde rotation with small obliquities, which is consistent with semi-collisional or 
collisional accretion. However, no firm conclusions can be drawn from such a small data set 
and we cannot rule out that the low obliquities of Earth and Mars are coincidental. 

Terrestrial planet formation in the asteroid belt was interrupted when growing planets 
became massive enough to gravitationally perturb the local population, causing bodies to 
collide with increased energy, ending accretion and commencing fragmentation. Evidence 
from Vesta's crust (e.g. Chapman 1986) and recent models of collisional evolution in the 
astroid belt (Gil-Hutton 1997; Bottke et al. 2005) suggest that the largest asteroids have 
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survived un-shattered and that they experienced very little collisional evolution. Their cur- 
rent spin properties may therefore still contain some information about their primordial spin 
state and hence clues about the formation of protoplanets (Davis et al. 1989; Bottke et al. 
2005). The two most massive asteroids, Ceres and Vesta, both exercise prograde rotation 
with periods of 9.1 and 5.3 hours respectively. Ceres' spin axis has a 12° inclination with 
respect to the normal of the ecliptic (Thomas et al. 2005) and Vesta's spin axis inclination to 
the normal of the ecliptic is ~ 40° (Drummond et al. 1998). The spin properties of Ceres and 
Vesta might therefore be indicative of semi-collisional or collisional accretion in the asteroid 
belt. 

Kuiper Belt objects (KBOs) grew mainly by planetesimal accretion. The formation 
time for Pluto-sized KBOs is comparable to the time required for a collisional cascade to set 
in, grinding initially kilometer-sized planetesimals to meters in size. If indeed a collisional 
cascade started by the time the largest KBOs formed, semi-collisional accretion could have 
dominated their formation. This may explain the intriguingly rapid spin of 2003 EL 61 , 
whose rotation period is only ~ 4 hours (Rabinowitz et al. 2006). However, the retrograde 
rotations of Pluto and 2003 EL@i (Brown et al. 2005) (assuming that it spins in the same 
direction as it is orbited by its largest satellite) conflict with this and tentatively suggest 
that semi-collisional accretion did not dominate their formation. 
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Fig. 1. — Collisionless and semi-collisional accretion. The protoplanet is represented by the 
filled black circle and its Hill radius is given by the solid black line. In the case of collisionless 
accretion (light grey) only planetesimals with impact parameters that allow direct collision 
with the protoplanet are accreted. In the semi-collisional and collisional regimes (dark grey) 
unbound planetesimals collide inside the Hill sphere of the protoplanet producing bound 
planetesimals which form a prograde accretion disk around the protoplanet. This enables 
the protoplanet to effectively accrete at its Hill radius. 
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Fig. 2. — Mean specific angular momentum < l z > in units of VlR 2 H accreted from a cold disk 
of planetesimals vs accretion radius R acc in Hill radii. The crosses indicate the results from 
our numerical integration and the dashed line corresponds to the limit in which the gravity 
of the protoplanet can be neglected, i.e. R acc >> Rh- The solid line shows the analytic 
solution valid for R acc « R H (Dones & Tremaine 1993a). 



